Thermoelectric and thermal transport in bilayer graphene systems 
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We numerically study the disorder effect on the thermoelectric and thermal transport for bilayer 
graphene under a strong perpendicular magnetic field. In the unbiased case, we find that the 
thermoelectric transport has similar properties as in the monolayer graphene, i.e., the Nernst signal 
has a peak at the central Landau level (LL) with the value of the order of fcs/e and changes sign near 
other LLs while the thermopower has an opposite behavior. We attribute this to the coexistence 
of particle and hole LLs around the Dirac point. When a finite interlayer bias is applied and a 
band gap is opened, it is found that the transport properties are consistent with those of a band 
insulator. We further study the thermal transport from electronic origins and verify the validity of 
the generalized Weidemann-Franz law. 

PACS numbers: 72.80.Vp; 72.10.-d; 73.50.Lw, 73.43.Cd 



I. INTRODUCTION 

Thermoelectric transport properties of graphene have 
attracted much recent experimental ItlliSi] and theoret- 
ical [l-Q attention. The thermopower (the longitudi- 
nal thermoelectric response) and the Nernst signal (the 
transverse response) in the presence of a strong magnetic 
field are found to be large, reaching the order of the quan- 
tum limit ks/e, where fc^ and e are the Boltzmann con- 
stant and the electron charge, respectively 1-3]. This 
has been attributed to the semi-metal type dispersion of 
graphene and/or in the vicinity of a quantum Hall liq- 
uid to insulator transition where the imbalance between 
the particle and hole types of carriers is significant. The 
thermoelectric effects are very sensitive to such an imbal- 
ance and become large in comparison with conventional 
metals. 

In our previous study on graphene in the presence 
of disorder and an external magnetic field @, we have 
shown that its thermoelectric transport properties are 
determined by the interplay of the unique band struc- 
ture, the disorder-induced scattering, the Landau quan- 
tization and the temperature. While the band structure 
and the magnetic field determine the Landau level (LL) 
spectrum, the broadening of each LL is controlled by the 
competition between disorder-induced scattering and the 
thermal activation. We find that all transport coeffi- 
cients are universal functions of Wl/Ef and kBT/Ep 
when both Wl and fc^T are much smaller than the Lan- 
dau quantization energy huc- Here Wl, Ep and T are 
the disorder-induced LL broadening, the Fermi energy 
and the temperature, respectively. When ksT <^ Wl, 
the thermoelectric conductivities vary as the density of 
states (and the particle-hole symmetry) is tuned by Ep 
from the center of the LL to the mobility gap. When 
ksT 3> Wl, thermal activation dominates and certain 
peak values for the thermopower Sxx or the Nernst signal 



Sxy reach universal numbers independent of the magnetic 
field or the temperature. While both Sxx and Sxy near 
high LLs {i/ 0) have similar behaviors as those in two- 
dimensional (2D) semiconductor systems displaying the 
integer quantum Hall effect (IQHE) p^| - [l3j . they rather 
have opposite behaviors around the central LL. Sxy has 
a peak while Sxx vanishes and changes sign at the Dirac 
point {Ep = 0). We have further argued that the unique 
behavior at the central LL is due to the coexistence of 
particle and hole LLs. As protected by the particle-hole 
symmetry, the contributions from particle and hole LLs 
cancel with each other exactly in the thermopower but 
superpose in the Nernst signal. The results for such a 
tight-binding analysis ar_e in good agreement with exper- 
imental observations [u- 



In this work, we extend our study to bilayer graphene 
which has two parallel graphene sheets stacked on top of 
each other as in 3D graphite (the AB or Bernal stacking) . 
While some common features are observed related to LLs 
with the same underlying particle-hole symmetry, bilayer 
graphene also demonstrates some interesting and differ- 
ent properties from monolayer graphene |14l4l9l | . The low 
energy dispersion of bilayer graphene can be effectively 
given by two hyperbolic bands Ck ~ ±A:^/(2to*) touching 
each other at the Dirac point {Ep — 0), i.e., the electrons 
or holes have a finite mass m* which is in contrast to the 
massless excitations in monolayer graphene. Another im- 
portant difference of bilayer graphene is the possibility to 
open up a band gap with a bias voltage, or a potential 
difference, applied between the two layers. This tunable 
gap system is advantageous to conventional semiconduc- 
tor materials, making bilayer graphene more appealing 
from the point of view of applications. The thermoelec- 
tric transport properties of bilayer graphene are also ex- 
pected to be interesting. The the thermopower of bi- 
layer graphene without a magnetic field has been consid- 
ered |20| . It is shown that as the density of states is also 
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of the pseudogap type without a biased voltage, one ex- 
pects that the relation for the thermopower S^x ~ T/Ep 
continues to hold. In addition, it is found that the room- 
temperature thermopower with a bias voltage can be en- 
hanced by a factor of 4 than monolayer graphene or un- 
biased bilayer graphene [lO], making it a more promis- 
ing candidate for future thermoelectric applications. Our 
study is to consider the thermopower and the Nernst ef- 
fect under a magnetic field. 

When an external magnetic field B is applied, as in 
graphene and other IQHE systems, electron states of 
bilayer graphene are quantized into Landau levels. As 
the band dispersion changes, these LLs follow a differ- 
ent quantization sequence En = zty''n{n — l)hajc with 
(jjc ^ B rather than ^/B for graphene. This has been 
confirmed by the theoretical \2l] and experimental [l^l 
studies on the quantum Hall effects, and further veri- 
fied by our numerical calculation [23|. Compared with 
graphene, though the massive nature of particles and hy- 
perbolic dispersion are different, the existence of the cen- 
tral LL (v = 0) and the associated chiral and particle- 
hole symmetries are preserved. Therefore, the study on 
the thermoelectric transport in bilayer graphene not only 
provides theoretical predictions for their properties, in 
particular, their dependence on disorder and magnetic 
field for this system, but also helps to verify our argu- 
ment on the central LL that its unique behavior is due to 
the chiral and particle-hole symmetries associated with 
the Dirac point. 

For such purposes, we carry out a numerical study of 
the thermoelectric transport in both unbiased and bi- 
ased bilayer graphene. We focus on studying the effects 
of disorder and thermal activation on the broadening 
of LLs and the corresponding thermoelectric transport 
properties. In the unbiased case, we indeed observe sim- 
ilar behaviors as in monolayer graphene for the central 
LL. Both the longitudinal and the transverse thermo- 
electric conductivities are universal functions of Wl/Ef 
and ksT/Ep and display different asymptotic behaviors 
in different temperature regions. The calculated Nernst 
signal has a peak at the central LL with heights of the 
order of fc^/e, and changes sign near other LLs, while the 
thermopower has an opposite behavior. A higher peak 
value is obtained comparing to graphene due to the dou- 
bled degeneracy. This confirms our argument that as the 
particle and hole LLs coexist only in the central LL, the 
thermopower vanishes while the Nernst effect has a peak 
structure. As before, we verify the validity of the semi- 
classical Mott relation, which is shown to hold in a wide 
range of temperatures. When a bias is applied between 
the two graphene layers, the thermoelectric coefficients 
exhibit unique characteristics quite different from those 
of unbiased case. Around the Dirac point, the trans- 
verse thermoelectric conductivity exhibits a pronounced 
valley with Uxy = at low temperature, and the ther- 
mopower displays a very large peak. We show that these 
features are associated with a band insulator, due to the 
opening of a sizable gap between the valence and con- 



ductance bands in biased bilayer graphene. In addition, 
we have calculated the thermal transport properties of 
electrons for both unbiased and biased bilayer graphene 
systems. In the biased case, it is found that the trans- 
verse thermal conductivity displays a pronounced plateau 
with — 0, which is accompanied by a valley 
This provides additional evidence for the band insulator 
behaviors. We further compare the calculated thermal 
conductivities with those deduced from the Wiedemann- 
Franz law, to check the validity of this fundamental rela- 
tion in graphene systems. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian. In Sec. Ill and Sec. 
IV, numerical results based on exact diagonalization and 
thermoelectric transport calculations are presented for 
unbiased and biased systems, respectively. In Sec. V, nu- 
merical results for thermal transport are presented. The 
final section contains a summary. 



II. MODEL AND METHODS 



We consider a bilayer graphene sample consisting of 
two coupled hexagonal lattices including inequivalent 
sublattices A, B on the bottom layer and A, B on the top 
layer. The two layers are arranged in the AB (Bernal) 
stacking [l^, HI] , where B atoms are located directly be- 
low A atoms, and A atoms are the centers of the hexagons 
in the other layer. Here, the in-plane nearest-neighbor 
hopping integral between A and B atoms or between A 
and B atoms is denoted by jab — IXb ~ ^o- For the 
interlayer coupling, we take into account the largest hop- 
ping integral between B atom and the nearest A atom 
= 7i, and the smaller hopping integral between an 

A atom and three nearest B atoms 7^^^ = 73. The values 
of these hopping integrals are taken to be 70 = 3.16 eV, 
71 = 0.39 eV, and 73 = 0.315 eV, as same as in Ref. [23| . 

We assume that each monolayer graphene has to- 
tally L„ zigzag chains with atomic sites on each 
chain [2g. The size of the sample will be denoted as 
N = Lx X Ly X Lz, where = 2 is the number of 
monolayer graphene planes along the z direction. We 
model charged impurities in substrate, randomly located 
in a plane at a distance d from the graphene sheet with 
long-range Coulomb scattering potentials 27]. This type 
of disorder is known to give more satisfactory results for 
transport prop erties of graphene in the absence of a mag- 
netic field [2^1 . When a magnetic field is applied perpen- 
dicular to the bilayer graphene plane, the Hamiltonian 
can be written in the tight-binding form 
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{ij)a {ij)a 



/i.e. 



(1) 



where cj^(c^^^) are creating operators on 



and cl„ 



A and B sublattices in the bottom layer. 
c!„(c^ fv) are creating operators on A and B sublattices 
in the top layer, tr is a spin index. The sum 
denotes the intralayer nearest-neighbor hopping in both 
layers, X)(ij)i(T stands for interlayer hopping between the 
B sublattice in the bottom layer and the A sublattice in 
the top layer, and X](ij)3(T stands for the interlayer hop- 
ping between the A sublattice in the bottom layer and 
the B sublattice in the top layer, as described above. 
The magnetic flux per hexagon (j) — a^j = ^ is pro- 
portional to the strength of the applied magnetic field 
B, where M is assumed to be an integer and the lat- 
tice constant is taken to be unity. For charged impu- 
rities, Wi — —^^J2a l/V ~ R-a)^ + d"^, where Ze is 
the charge carried by an impurity, e is the effective back- 
ground lattice dielectric constant, and and are the 
planar positions of site i and impurity a, respectively. All 
the properties of the substrate (or vacuum in the case of 
suspended graphene) can be absorbed into a dimension- 
less parameter = Z / {ehv f) , where vf is the Fermi 
velocity of the electrons. For simplicity, in the following 
calculation, we fix the values of distance d = Inm and 
impurity density as 1% of the total sites, and tune to 
control the impurity scattering strength. 

For the biased system, the two graphene layers gain 
different electrostatic potentials, and the corresponding 
energy difference is given by Ag — ^2 ~ where ei — 
— iA„, and ea = ^A„. The Hamiltonian can be written 



as: 



H = Hn 



J2'^i{CiaCza + e2cl„Cia). For illustrative 

i(T 

purpose, a relatively large asymmetric gap A ^ — O.I70 is 
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assumed, which is experimentally achievable 

In the linear response regime, the charge current in 
response to an electric field or a temperature gradient 
can be written as J = (tE -|- d(— VT), where a and d are 
the electrical and thermoelectric conductivity tensors, re- 
spectively. These transport coefficients can be calculated 
by Kubo formula once we obtain all the eigenstates of the 
Hamiltonian (in our calculation, axx is obtained based on 
the calculation of the Thouless number [1^). In practice, 
we can first calculate the T = conductivities aji{EF), 
and then use the relation 



aj^{EF,T) = J d<iaj^{<i){t- Ef)(^ 



dm 

de 



to obtain the finite temperature electrical and ther- 
moelectric conductivity tensors. Here, f{x) = 
\l[g(^-EF)/kBT ^ |-]^g Fermi distribution function. 
At low temperatures, the second equation can be ap- 
proximated as 



aji{EF,T) 



3e 



de 



(3) 



<i=Ef 



which is the semiclassical Mott relation [T^, The 
thermopower and Nernst signal can be calculated subse- 
quently from (29} 
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Pyx^yx ■) 
Pyx^xx ■ 



(4) 



The thermal conductivity, measuring the magnitude 
of the thermal currents in response to an applied tem- 
perature gradient, includes electron and phonon contri- 
butions. In our numerical calculations, phonon-derived 
thermal conductivity is omitted. The electronic ther- 
mal conductivities Kji at finite temperature assume the 
forms 11311 



J dea,,{e){e-EFf 



df{e) 



Taj,{EF.T)<jJ^\EF. T)aj^{EF, T). (5) 



For diffusive electronic transport in metals, it is well 
known that the Wiedemann-Franz law is satisfied be- 
tween the electrical conductivity a and the thermal con- 
ductivity K of electrons (30j : 



(6) 



where L is the Lorentz number and takes a constant 
value: L=^{^f. 



III. THERMOELECTRIC TRANSPORT IN 
UNBIASED BILAYER GRAPHENE 

We first show calculated thermoelectric conductivities 
at finite temperatures for unbiased bilayer graphene. As 
seen from FiglIJa) and (b), the transverse thermoelectric 
conductivity axy displays a series of peaks, while the lon- 
gitudinal thermoelectric conductivity Uxx oscillates and 
changes sign at the center of each LL. At low temper- 
atures, the peak of axy at the central LL is higher and 
narrower than others, which indicates that the impurity 
scattering has less effect on the central LL. These re- 
sults are qualitatively similar to those found in mono- 
layer graphene [1] due to the similar particle-hole sym- 
metry in both cases, but some obvious differences exist. 
Firstly, the peak values of axy at the central LL is larger 
than that of monolayer graphene. Secondly, at low tem- 
perature, axy splits around Ep — ±0.4670, which can 
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FIG. 1: (color online). Thermoelectric conductivities at fi- 
nite temperatures of bilayer graphene. (a)-(b) axy{EF,T) 
and axx{EF,T) as functions of the Fermi energy at differ- 
ent temperatures, (c) shows the temperature dependence of 
oixyiEp ,T) for monolayer and bilayer graphene. (d) com- 
pares the results from numerical calculations and from the 
generalized Mott relation at two characteristic temperatures, 
UbT/Wl = 0.05 and UbT/Wl = 1. The system size is taken 
to be = 96 X 48 X 2, magnetic flux (f> = 27r/48, and disorder 
strength = 0.3 (we consider uniformly distributed positive 
and negative charged impurities within this strength) with 
WlHo = 0.0376. 



be understood as due to the presence oi v = ±8 Hall 
plateau by lifting subband degeneracy. In Figdjc), we 
find that axy shows different behavior depending on the 
relative strength of temperature fc^T and the width of 
the central LL Wl{Wl is determined by the full- width at 
half- maximum of the a^x peak). When kgT Wl and 
Ep <^ Wl, axy shows linear temperature dependence, 
indicating that there is a small energy range where ex- 
tended states dominate, and transport fall into the semi- 
classical Drude-Zener regime. When Ep is shifted away 
from the Dirac point, the low temperature electron exci- 
tation is gapped related to Anderson-localization. When 
ksT becomes comparable to or greater than Wl, the axy 
for all LLs saturates to a constant value ^.bAkse/h. This 
matches exactly the universal number {\n2)kBe/h pre- 
dicted for the conventional IQHE systems in the case 
where thermal activation dominates [l^ [l3l |. with an 
additional degeneracy factor 8. The saturated value of 
Uxy in bilayer graphene is exactly twice than that of the 
monolayer graphene, as shown in Figlljc) in accordance 
with the eightfold degeneracy from valley, spin and layer 
degree of freedoms [2l|, [22] . 

To examine the validity of the semiclassical Mott rela- 
tion, we compare the above results with those calculated 
from Eq.(l3]), as shown in Fig[T{d). The Mott relation is 
a low-temperature approximation and predicts that the 
thermoelectric conductivities have linear temperature de- 
pendence. This is in agreement with our low-temperature 



FIG. 2: (color online). Thermoelectric conductivities at finite 
temperatures of biased bilayer graphene. (a)-(b) axyi^Ep ,T) 
and axxiEp ,T) as functions of the Fermi energy at different 
temperatures, (c) Compares the results from numerical cal- 
culations and from the generalized Mott relation at two char- 
acteristic temperatures, ftsT/Ac, = 0.05 and ksT j — 0.2. 
Here asymmetric gap Ag = O.I70. The system size is taken 
to be A'' = 96 X 48 X 2, magnetic flux (j> = 27r/48, and disorder 
strength = 0.3. 



results, which proves that the semiclassical Mott relation 
is asymptotically valid in Landau-quantized systems, as 
suggested in Ref. [l2j |. 
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(b) bilayer graphene 
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(d) biased bilayer graphene 
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FIG. 3: (color online). The thermopower Sxx and the Nernst 
signal Sxy as functions of the Fermi energy in (a)-(b)bilayer 
graphene, (c)-(d)biased bilayer graphene at different temper- 
atures. All parameters in this two systems are chosen to be 
the same as in Fig. 1 and Fig. 2, respectively. 
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IV. THERMOELECTRIC TRANSPORT IN 
BIASED BILAYER GRAPHENE 

For biased bilayer graphene, we show results of axx and 
a^y at finite temperatures in Fig. [5] Here we see that a^y 
demonstrates a pronounced valley, in striking contrast to 
the unbiased case with a peak at the particle-hole sym- 
metric point Ef = 0. This behavior can be understood 
as due to the split of the valley degeneracy in the cen- 
tral LL by an opposite voltage bias added to the two 
layers. This is consistent with the opening of a sizable 
gap between the valence and conduction bands. More 
oscillations are observed in the higher LLs comparing to 
the unbiased case, in consistent with the further lifting 
of the LL degeneracy, axx oscillates and changes sign 
around the center of each split LL. In Figl2ljc), we also 
compare the above results with those calculated from the 
semiclassical Mott relation using Eq.Q. Here the Mott 
relation is shown to remain valid at low temperature. 

We further calculate the thermopower Sxx and the 
Nernst signal Sxy using Eq. which can be directly 
determined in experiments by measuring the responsive 
electric fields. In Fig. [3lja)-(b), we show results of Sxx 
and Sxy in unbiased bilayer graphene. As we can see, 
Sxy (Sxx) has a peak at the central LL (the other LLs), 
and changes sign near the other LLs (the central LL), 
similar to the case of monolayer graphene foj. This oscil- 
latory feature has been observed experimentally [31]. In 
our calculation, the peak value of Sxx at n = — 1 LL is 
found to be U^iV/K (note that ke/e = 86.17^1^//^ ) for 
ksT = 0.05Wl and 26fiV/K for keT = O.IW^L: which is 
in good agreement with the measured value [3l| . At zero 
energy, both pxy and axx vanish, leading to a vanishing 
Sxx- Around the zero energy, because PxxCtxx and PxyOxy 
have opposite signs, depending on their relative magni- 
tudes, Sxx could either increases or decreases when Ep 
is increased passing the Dirac point. In bilayer graphene, 
we find that Sxx is always dominated by PxyfXxy, conse- 
quently, Sxx decreases to negative value as Ep passing 
zero. We find that the peak value of Sxx in the cen- 
tral LL is ±6fiV/K at ksT = 0.05VFl. On the other 
hand, Sxy has a peak structure at zero energy, which 
is dominated by Pxx<^xy The peak value is 'i2pV/K at 
ksT = 0.05Wl- These results are in good agreement 
with the experiments. 

In Fig. [31[c)-(d), we show the calculated Sxx and Sxy 
in biased bilayer graphene system. As we can see, Sxy 
(Sxx) has a peak around zero energy (the other LLs), and 
changes sign near the other LLs (zero energy). In our cal- 
culation, Sxx is dominated by PxxO^xx^ which is different 
from the unbiased bilayer graphene. At low temperature, 
the peak value of Sxx around zero energy keeps almost 
unchanged around ±l81pV/ K , which is much larger than 
that of unbiased case. With the increase of temperature, 
the peak height increases to ±396pV/K at fc^T — 0.5Ag. 
Theoretical study [l^] indicates that, the large magni- 
tude of Sxx is mainly a result of the energy gap. On the 
other hand, Sxy has a peak structure around zero energy. 
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FIG. 4: (color online) 
:,y{EF,T) and 



(a)-(b) Thermal conductivities 
iixx{EF,T) as functions of the Fermi en- 
ergy in bilayer graphene at different temperatures, (c)- 
(d) Compares the thermal conductivity as functions of the 
Fermi energy from numerical calculations and from the 
Wiedemann-Franz Law at two characteristic temperatures. 
The parameters chosen here are the same as in Fig. 1. 



which is dominated by OxyPxx- With axx ^ 2e^//i near 
Ep = 0, we find that the peak height is 198pV/K at 
ksT = O.lAg, which is larger than that of unbiased case. 



(a) biased bilayer graphene 
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FIG. 5: (color online). (a)-(b) Thermal conductivities 
Hxy{EF,T) and HxxiEpjT) as functions of the Fermi en- 
ergy in biased bilayer graphene at different temperatures, 
(c)-(d)Compares the thermal conductivity as functions of 
the Fermi energy from numerical calculations and from the 
Wiedemann-Franz Law at two characteristic temperatures. 
The parameters chosen here are the same as in Fig. 2. 
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V. THERMAL CONDUCTIVITY FOR 
UNBIASED AND BIASED BILAYER 
GRAPHENE SYSTEMS 

We now focus on thermal conductivities. In Fig. |4l 
we show results of the transverse thermal conductivity 
Kxy and the longitudinal thermal conductivity Kxx for 
unbiased bilayer graphene at different temperatures. As 
seen from Figllja) and (b) , K^y exhibits two flat plateaus 
away from the center of the central LL. At low tem- 
peratures, the transition between these two plateaus is 
smooth and monotonic, while at higher temperatures, 
Kxy exhibits an oscillatory feature at ksT = 0.5Wl be- 
tween two plateaus. On the other hand, k^x displays a 
peak near the center of the central LL, while its peak 
value increases quickly with T. To test the validity of 
the Wiedemann- Franz Law, we compare the above results 
with ones calculated from Eq.® as shown in FigHJc) and 
(d). The Wiedemann- Franz Law predicts that the ratio 
of the thermal conductivity k to the electrical conductiv- 
ity (T of a metal is proportional to the temperature. This 
is in agreement with our low-temperature results, while 
deviation is seen at higher T. 

In Fig. [5l we show the calculated thermal conduc- 
tivities Kxx and Kxy for biased bilayer graphene. As 
seen from FiglSfa) and (b), around the zero energy, a 
flat region with Kxy = is found at low temperatures, 
which is accompanied by a valley in Kxx- These features 
are clearly in contrast to those of unbiased case due to 
the asymmetric gap between the valence and conduction 
bands. When temperature increases to ksT — 0.2Ag, 
the plateau with Kxy — disappears, while Kxx displays 
a large peak. In FiglSfc) and (d), we also compare above 
results with those calculated from the Wiedemann-Franz 
Law using Eq.®. Due to the presence of energy gap, we 
find that the Wiedemann-Franz Law is not valid in the 
biased bilayer graphene. 

VI. SUMMARY 

In summary, we have numerically investigated the 
thermoelectric and thermal transport in unbiased bi- 
layer graphene based on the tight-binding model in the 
presence of both disorder and magnetic field. We find 
that the thermoelectric conductivities display different 
asymptotic behaviors depending on the ratio between 
the temperature and the width of the disorder-broadened 
Landau levels (LLs) , similar to those found in monolayer 
graphene. In the high temperature regime, the transverse 
thermoelectric conductivity axy saturates to a universal 
quantum value 5.54fcse//i at the center of each LL, and 



it has a linear temperature dependence at low tempera- 
tures. The calculated Nernst signal Sxy has a peak at the 
central LL with heights of the order of fcs/e, and changes 
sign at the other LLs, while the thermopower Sxx has an 
opposite behavior. These results are in good agreement 
with the experimental observation[3l|. The validity of 
the semiclassical Mott relation between the thermoelec- 
tric and electrical transport coefficients is verified in a 
range of temperatures. The calculated transverse ther- 
mal conductivity Kxy exhibits two plateaus away from the 
band center. The transition between this two plateaus is 
continuous, which is accompanied by a pronounced peak 
in longitudinal thermal conductivity Kxx- The validity of 
the Wiedemann-Franz Law between the thermal conduc- 
tivity K and the electrical conductivity a is only verified 
at very low temperatures. 

We further discuss the thermoelectric transport of bi- 
ased bilayer graphene. When a bias is applied to the two 
graphene layers, the thermoelectric coefficients exhibit 
unique characteristics different from those of unbiased 
case. Around the Dirac point, transverse thermoelectric 
conductivity exhibits a pronounced valley with axy — 
at low temperatures, and the thermopower displays a 
large magnitude peak. Furthermore, the transverse ther- 
mal conductivity has a pronounced plateau with Kxy = 0, 
which is accompanied by a valley These are in 

consistent with the opening of sizable gap between the va- 
lence and conductance bands in biased bilayer graphene. 

We mention that in our numerical calculations, the 
magnetic field is much stronger than the ones one can re- 
alize in the experimental situation, as limited by current 
computational capability. However, the asymptotic be- 
haviors we obtained is robust and applicable to weak field 
limit since it is determined by the topological property 
of the energy band as clearly established for monolayer 
graphene [0]. 

Acknowledgments 

This work is supported by the DOE Office of Basic 
Energy Sciences under grant DE-FG02-06ER46305 (RM, 
DNS), and the U.S. DOE through the LDRD program 
at LANL (LZ), the NSF Grant DMR-0906816 (RM). 
We also thank partial support from Princeton MRSEC 
Grant DMR-0819860, the NSF instrument grant DMR- 
0958596 (DNS), the NSFC Grant No. 10874066, the 
National Basic Research Program of China under Grant 
Nos. 2007CB925104 and 2009CB929504 (LS), and the 
doctoral foundation of Chinese Universities under Grant 
No. 20060286044 (ML). 



[1] Y.M. Zuev, W. Chang, and P. Kim, Phys. Rev. Lett. [2] P. Wei, W. Bao, Y. Pu, C. N. Lau, and J. Shi, Phys. Rev. 
102, 096807 (2009). Lett. 102, 166808 (2009). 



7 



[3] J.G. Checkclsky and N.P. Ong, Phys. Rev. B 80, 

081413(R) (2009). 
[4] A. H. Castro Ncto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

[5] E. H. Hwang, E. Rossi, and S. Das Sarma, Phys. Rev. B 

80, 235415 (2009). 
[6] T. Lofwander and M. Fogelstrom, Phys.Rev. B 76, 

193401 (2007). 

[7] B. Dora and P. Thalmeier, Phys. Rev. B 76, 035402 
(2007). 

[8] X.-Z. Yan, Y. Romiah, and C. S. Ting, Phys. Rev. B 80, 
165423 (2009). 

[9] L. Zhu, R. Ma, L. Sheng, M. Liu, and D. N. Sheng, Phys. 

Rev. Lett. 104, 076804 (2010). 
[10] S. M. Girvin and M. Jonson, J. Phys. C 15, L1147(1982). 
[11] P. Stfeda, J. Phys. C 16, L369 (1983). 
[12] M. Jonson and S.M. Girvin, Phys. Rev. B 29, 1939 

(1984). 

[13] H. Oji, ,J. Phys. C 17, 3059 (1984). 

[14] E. McCann, Phys. Rev. B 74, 161403(R) (2006). 

[15] H. Min, B. Sahu, S. K. Banerjee, and A. H. MacDonald, 

Phys. Rev. B 75, 155115 (2007). 
[16] E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 

Peres, J. M. B. Lopes dos Santos, J. Nilsson, F. Guinea, 

A. K. Geim, and A. H. Castro Neto, Phys. Rev. Lett. 99, 

216802 (2007). 

[17] J. B. Oostinga, H. B. Hcersche, X. Liu, A. F. Morpurgo, 
and L. M. K. Vandersypen, Nat. Mater. 7, 151 (2008). 

[18] Y. Zhang, T.-T. Tang, C. Girit, Z. Hao, M. C. Martin, A. 
Zettl, M. F. Crommie, Y. R. Shen, and F. Wang, Nature 



459, 820 (2009). 
[19] K. F. Mak, C. H. Lui, J. Shan, and T. F. Heinz, Phys. 

Rev. Lett. 102, 256405 (2009). 
[20] L. Hao and T. K. Lee, Phys. Rev. B 81, 165445 (2010). 
[21] E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 

(2006). 

[22] K. S. Novoselov, E. McCann, S. V. Morozov, V. L Fal'ko, 

M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin and A. 

K. Geim, Nature Phys. 2, 177 (2006). 
[23] R. Ma, L. Sheng, R. Shen, M. Liu and D. N. Sheng, Phys. 

Rev. B 80, 205101 (2009); R. Ma, L. Zhu, L. Sheng, M. 

Liu, D. N. Sheng, Europhys. Lett. 87, 17009 (2009). 
[24] S. B. Trickey, F. Miiller-Plathe, and G. H. F. Diercksen, 

Phys. Rev. B 45, 4460 (1992). 
[25] K. Yoshizawa, T. Kato, and T. Yamabe, J. Chem. Phys. 

105, 2099 (1996); T. Yumura and K. Yoshizawa, Chem. 

Phys. 279, 111 (2002). 
[26] D.N. Sheng, L. Sheng, and Z.Y. Weng, Phys. Rev. B 73, 

233406 (2006). 

[27] S. Adam and S. Das Sarma, Solid State Communications 

146,356 (2008). 
[28] S. Adam, E. H. Hwang, V. M. Galitski, and S. Das 

Sarma, Proc. Natl. Acad. Sci. USA 104, 18392 (2007). 
[29] Different literatures may have a sign difference due to 

different conventions. 
[30] J. M. Ziman, Electrons and Phonons (Oxford University 

Press, 

[31] S.G. Nam, D.K. Ki, H.J. Lee, Phys. Rev. B 82, 245416 
(2010). 



